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Abstract 



Adopting the hypothesis about the exact cancellation of vacuum condensates contributions to the 
ground state energy in particle physics to the leading order in graviton-mediated interactions, we 
argue that the observable cosmological constant can be dynamically induced by an uncompensated 
quantum gravity correction to them. To start with, we demonstrate a possible cancellation of the 
. quark-gluon condensate contribution to the total vacuum energy density of the Universe without 

taking into account the graviton-mediated effects. In order to incorporate the latter, we then 
calculate the leading-order quantum correction to the classical Einstein equations due to metric 
fluctuations induced by the non-perturbative vacuum fluctuations of the gluon and quark flelds in 
the quasiclassical approximation. It has been demonstrated that such a correction to the vacuum 
energy density has a form ea ~ GAq(^j-,, where G is the gravitational constant, and Aqcd is the 
QCD scale parameter. We analyze capabilities of this approach based on the synthesis between 
quantum gravity in quasiclassical approximation and theory of non-perturbative QCD vacuum for 
quantitative explanation of the observed Dark Energy density. 
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I. INTRODUCTION 



The existence of mysterious Dark Energy which drives current accelerated expansion of 
the Universe is confirmed in many cosmological observations so far, e.g. in studies of the 
type la Supernovae [1], cosmic microwave background anisotropies [2|, large scale structure 
[3| etc. The hypothesis about the time-independent cosmological constant called A-term is 
an essential part of the Standard Cosmological Model known as the Cold Dark Matter with 
A-term Model (or ACDM) and agrees well with all observational data collected so far at 
the current level of experimental accuracy. However, theoretical origin of the cosmological 
constant has not been properly understood yet. 

Nowadays, the Dark Energy problem remains one of the major unsolved problem of 
theoretical physics 01 . On the way of searching for possible solutions of this problem many 
various pathways were explored during last few decades referring to e.g. new exotic forms 
of matter (e.g. "quintessence" jsj, "phantom" jof etc), holographic models 0], string theory 
landscape js], lof, Born-Infeld quantum condensate 10 1, modified gravity approaches 11, 12 



etc. For a comprehensive overview of existing theoretical models and interpretations of the 
Dark Energy, see e.g. Refs. [l3-16[ and references therein. In this work we are primarily 
focused on the class of those approaches which are based upon conventional quantum field 
theory and standard quantum gravity in quasiclassical approximation. 

The traditional identification of the A-term satisfying the equation of state Pa = — 
with vacuum pressure Pa and energy density suffers from a large gap of knowledge on 
non-perturbative dynamics of the ground state in particle physics. In particular, individual 
vacuum condensates e.g. those which are responsible for the chiral (quark-gluon condensate) 
and gauge (Higgs condensate) symmetries breaking in the Standard Model contribute to the 
vacuum energy of the Universe individually exceeding the observable value of the Dark 
Energy by many orders of magnitude in absolute value [13]. This situation (sometimes 
referred to as the "Vacuum Catastrophe" in the literature) requires extra hypotheses about 
(partial or complete) compensation of vacuum condensates of different types to the net 
vacuum energy density of the Universe (see e.g. Ref. jof). Note that most of the models 
developed in the literature adopt the same point of view trying to cancel or suppress short 
distance vacuum fluctuations in one way or another (for a review on such models, see e.g. 
Ref. [l3| and references therein). Still, a dynamical mechanism for such gross cancellations 
and corresponding major fine-tuning of vacuum parameters is not known at the moment 
and is a subject of intensive studies in the vast literature (see e.g. Ref. [11]). 

Historically, about 45 years ago Ya. B. Zeldovich has pointed out in Ref. [l^ that cos- 
mological constraints on the value of the A-term density to a good accuracy correspond to 
a simple purely mass dimensional estimate e\ ~ Gm^ (in natural units h = c = 1), where 
G = Mpf is the gravitational constant and m is some characteristic mass of light elementary 
particles (or light hadrons known at that time). For the first time, he proposed an interpreta- 
tion of the A-term as an effect of gravitational interactions of virtual particles in the physical 
vacuum. Almost at the same time, A. D. Sakharov noticed in Ref. jlof that, indeed, extra 
terms describing an effect of graviton exchanges between identical particles (e.g. bosons in 
the ground state) should appear in the right hand side of Einstein equations averaged over 
their quantum ensemble. Below we refer to this approach as the Zeldovich-Sakharov (ZS) 
scenario. 

One of the well-known and naive representations of the Zeldovich relation ~ Gm^ 
through the basic fundamental constants, the minimal (typical hadron scale) and maximal 
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(Planck scale) scales of fundamental particle physics, has been proposed by N. S. Kardashev 
in Ref. jllj. It has the following form 



where ~ 138 MeV is the pion mass |22|, Mpi = G'^^^ ~ 1.22 ■ 10^^ MeV is the 



Planck mass. It is worth to notice that the representation (11.11) turns out to be numerically 
very close to the most recent WMAP data, well within experimental error bars, e^^^ = 
(3.0±0.7) X 10"'^^ MeV^ [2*]. Such a remarkable numerical coincidence of the naive Zeldovich- 
Kardashev (ZK) formula (II. ip to the current cosmological observations seems to be almost 
too good to be just an accident [2^. This situation yet remains puzzling and requires a 
closer look into possible physical reasons for that. A deeper theoretical understanding of 
those reasons is the major goal of our current work. 

It is well-known that the pion mass is an object of essentially non-perturbative Quantum 
Chromodynamics (QCD) driven by properties of the quark-gluon vacuum condensate [13], 
the strongest non-perturbative vacuum sub-system at the minimal energy scales known in 
particle physics. So the numerical coincidence of the ZK result (11.11) to the observational 
data if not accidental poses a question about a possible role of strong non-perturbative 
vacuum fluctuations of quark and gluon fields in the A-term generation in the framework of 
ZS scenario described above (TH, 20 . 



Among existing quantum field theory approaches describing the observable A-term den- 
sity value, the most promising and successful ones are recently proposed by Klinkhamer and 
Volovik in Ref. j27j and by Urban and Zhitnitsky in Ref. jlsj . These approaches attempt to 
interpret the positive Dark Energy density as a small but non- vanishing effect of gravitating 
non-perturbative QCD vacuum fluctuations in a non-trivial background of expanding Uni- 
verse. The first approach jl^] is based upon a generic "g-theory" operating with a conserved 
microscopic q value, whose statics and dynamics are studied at macroscopic scales and which 
can, in principle, be identified with the gluon condensate in QCD. In this case it has been 
shown explicitly that the gravitating non-perturbative QCD vacuum fluctuations dynami- 
cally generate a nonzero limiting value of the vacuum energy density in the nonequilibrium 
context of the expanding Universe. The second approach [28| focuses on dynamics of the 
ghost fields in the low energy (chiral) QCD. In particular, it was shown that the Veneziano 
ghost, being unphysical in the usual Minkowski QFT, exhibits non-vanishing physical effect 
in the expanding universe leading to a positive vacuum energy density with a time- dependent 
equation of state. Both approaches arrive at similar order-or-magnitude estimates for the 
A-term density in the expanding Universe in terms of the maximal Mpi and minimal Aqcd 
scales of particle physics whose existence is required by the unitarity condition of an under- 
lined quantum theory [28|. Similar QCD-based ideas of the Dark Energy origin from the 



effective interacting gluon condensate were previously explored in Refs. [25[ . Also, quantum 
effects in the Born-Infeld fields condensation very similar to the gluon condensation in QCD 
(realized by means of the trace anomaly) which could potentially play an important role 
in the dynamical A-term generation were discussed in Ref. [lo|. Another interesting claim 
was made in Ref. jl^, where it was demonstrated that an extra contribution to the vacuum 
energy density of the form ~ GAqq-^ may also originate from a consideration of QCD in 
spacetime with torsion in the framework of the Einstein-Cartan theory with minimal fermion 
coupling to torsion. Such promising QCD-based interpretations of the Dark Energy give us 
a strong motivation for further investigations in this direction and requires development of 
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a rigorous self-consistent theoretical framework based upon quasiclassical (semiquantum) 
gravity and theory of non-perturbative QCD vacuum fluctuations, in their existing standard 
formulation. 

The basic ideas of this paper can be condensed into a few sentences as follows. In the 
zeroth order in metric fluctuations, the net energy density of non-perturbative QCD vacuum 
fluctuations to the energy density of the Universe can be, in principle, compensated via a 
macroscopic spatially-homogeneous modes of the cosmological Yang-Mills fields discussed 
by us recently in Ref. [29| or a microscopic dynamical mechanism of the QCD vacuum self- 
tuning. The latter possibility will be considered below in Section 11^. Further, in Section 
III and IV we demonstrate that the non-perturbative vacuum fluctuations of quark and 
gluon fields dynamically induce the fluctuations of the background metric in the expanding 
Universe through a coupling of gravity to the gluon field via the trace anomaly. Such a 
coupling leads to an extra correction term in the energy-momentum tensor which is linear 
in the graviton field (as the first-order quantum gravity correction in the quasiclassical ap- 
proximation). The existing (static) instanton theory of the non-perturbative QCD vacuum 
|3l| (for a review on the QCD theory of instantons see e.g. Ref. |i32j), however, does not 
enable us to constrain the precise numerical value of the QCD contribution to the A-term 
density fully dynamically, based on the First Principles only. Nevertheless, it is possible, 
in fact, to find some specific conditions and theoretical constraints under which the QCD 
vacuum theory leads to the observed A-term density value. Two of these conditions inherent 
to QCD, the presence of the conformal anomaly and a strong coupling theory, have already 
been emphasized e.g. in Refs. |23|, l33|. Besides, as will be shown in Section V and VI, under 
certain phenomenological constraints known from the chiral QCD theory and lattice simu- 
lations the extra quantum gravity correction term indeed gives rise to an additional positive 
non-vanishing A-term-type contribution to the vacuum energy density of the Universe close 
to the observable value, in accordance with the Sakharov hypothesis [201]. 



II. QCD VACUUM ENERGY: ZEROTH ORDER IN METRIC PERTURBATIONS 

The ground state in QCD is typically characterized by non-vanishing condensates of 
strongly interacting quarks and gluons commonly referred to as the quark- gluon condensate 
responsible for the confined phase of quark matter. In the framework of the instanton liq- 
uid models, the topological (or instanton) modes of the quark-gluon condensate are given 
essentially by the strong non-perturbative fluctuations of the gluon and sea (mostly, light) 
quark fields which are induced in processes of quantum tunneling of the gluon vacuum 
between topologically different classical states |32|. The topological instanton-type contri- 
bution e^acitop) to the energy density of the QCD vacuum is one of its main characteristics 



24j | and can be derived from the well-known trace anomaly relation [34 



^ Possible mechanisms for compensation of other perturbative contributions to the net vacuum energy density 
of the Universe from virtual fermions and bosons (e.g. Higgs condensates, graviton condensate etc) 



typically refer to Supersymmetry and high-scale Grand Unified Theories [l3|, |30| and will not be discussed 
further in this paper. 
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where T^^J-qcd) trace of the energy-momentum tensor, = Qg/An, iriq are the hght 

current sea quark masses, and F^^ is the gluon field strength tensor in the standard nor- 
malisation. A vacuum average of the trace (0|T?^qcj-))|0) = 4e:„ac(to£) finally leads to the 



well-known formula for the instanton energy density (see e.g. Refs. |35l. l36l|) 
;|(0|:^i^»(a;)Ff(x):|0) + i 

OZ TT 4 



vac[top) = -7^(0| : (x) : |0) + ^ (0| : m^uu : |0) + (0| : m^dd : |0) 



+ (0| : nisss : |0) ~ -(5 ± 1) x 10^ MeV^ , (2.2) 

composed of gluon and light sea m, d, s quark contributions. This is the saturated (maximal) 
value of the topological contribution to the QCD vacuum energy density, while possible 
refinements to it are dependent on poorl y k nown non-perturbative long-range Yang-Mills 
dynamics and were discussed e.g. in Ref. |36|. 

Clearly, contributions of a different physical nature should compensate the QCD contri- 
bution f l2.2p to the vacuum energy of the Universe since its value by far is not compatible 
with the cosmological observations and data on the A-term density value [2,]. On the other 
hand, within the general problem of vacuum condensates cancellation and corresponding 
fine tuning of vacuum substructures, the contribution (12. 2p has a special status. Various 
existing cancellation mechanisms refer essentiall y to an unknown high-scale physics beyond 



the Standard Model e.g. to Supersymmetry [13|, |30[. However, they cannot be applied for a 
compensation of the specifically non-perturbative and low-energy QCD contribution given 
by Eq. (12. 2p . This issue forces us to look for alternative cancellation mechanisms. 

One of the possible ways to eliminate the microscopic QCD vacuum contribution fl2.2p 
to the vacuum energy density of the Universe was discussed earlier by us in Ref. [29] in- 
troducing the hypothesis about the existence of the cosmological macroscopic Yang-Mills 
fields in early Universe. In particular, it was claimed that the negative QCD contribu- 
tion (12. 2 p can, in principle, be canceled by a positive constant contribution generated by 
spatially-homogeneous modes of such a field. Corresponding quasiclassical solution for these 
modes is exact, it necessarily takes into account interactions with the QCD vacuum (vacuum 
polarisation effects) and gives rise to the spatially-homogeneous finite-time instantons. 

In this paper, we consider another approach to the compensation of the topological 
QCD contribution assuming that there might be extra contributions from the spatially- 
inhomogeneous modes of a different nature to the QCD vacuum density at the QCD energy 
scale ~ Aqcd besides the topological (instanton) ones given by Eq. (12. 2 p and spatially- 



homogeneous ones predicted in Ref. 29|. In fact, the importance of such extra long-range 
quantum-wave fluctuations of the hadronic vacuum has been emphasized earlier, e.g. in 
Ref. [37|. Within the alternative microscopic approach to compensation of the QCD con- 
tribution, the total QCD vacuum energy density may turn into zero due to a fine-tuning 
of the QCD vacuum parameters corresponding to quantum-topological and quantum-wave 
fluctuations in a vicinity of the Aqcd scale. Let us explore such a possibility in detail. 
The quantum-topological fluctuations contributing to the QCD vacuum energy density 

12. 2p exist at typical space-time scales ~ Ig which satisfy the following approximate inequality 

32| 

lg{ra^n) ^ (1500 MeV)-\ ~ (500 McV)"^ . 

Here, the values of Ig(min), ig{max), which can be estimated e.g. in the lattice QCD framework, 
are interpreted as the minimal and maximal length scales of the non-perturbative gluon fleld 
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fluctuations and approximately correspond to boundaries of the light resonances region in 
the hadron spectrum [2^^. 

In general, besides topological (instanton-type) fluctuations contributing to Eq. fl2.2p . 
there are other two types of quantum-wave QCD vacuum fluctuations: (1) perturbative 
fluctuations of gluon and quark fields at smaller length scales / < lg(min), and (2) vacuum 
fluctuations corresponding to collective wave motion of gluon and quark fields at the same 
scales as the instanton ones fl2.3p [37]. As was mentioned earlier, the problem of compen- 
sation of the short-range perturbative QCD fluctuations, along with all other high energy 
vacuum subsystems corresponding to e.g. zero-point fluctuations of fundamental fields and 
Higgs-type condensates, is the subject of a supersymmetric "Theory of Everything" and 
therefore is not discussed here. Meanwhile, quantum-wave fluctuations of the second type 
(from now on, we denote them as the collective ones) have quantum numbers of hadrons 
with masses rrih < l~} . v Their renormalized Lorentz-invariant contribution to the net 
QCD vacuum energy density is then expressed in terms of the light hadron masses and the 
universal cut-off parameter /i ~ ^^(^m) (P^^yi^S ^ ^^^^ renormalisation scale) as follows 



where Jb and Jm {friB and mA/) are the spins (masses) of respective baryon and meson 
degrees of freedom, respectively. 

Note, both topological and collective fluctuations have a non-perturbative nature. Strictly 
speaking, the expression (12. 2 p based upon the trace anomaly relation (12. ip should contain 
both the quantum-topological and collective contributions. Their separate consideration, 
however, is necessary since they have a completely different structure [37| while there is no a 
self-consistent theory of the non-perturbative QCD vacuum dynamics which could enable us 
to extract the collective- wave contribution (12. 4p from the general expression (12. 2p . To this 
end, we formally remove the quantum-wave fluctuations from Eq. (12. 2 p taking the saturated 
(maximal) value for the topological contribution as was formally denoted by the symbols of 
normal ordering. 

The most important observation here is that taking into account only metastable hadronic 
degrees of freedom in Eq. (12. 4p - the baryon octet B = {N, A, S, S} and the pseudoscalar 
nonet M = {vr, K, r], r]'} - we obtain a meaningful result, namely, eyac{top) + £vac(h) = for 
a reasonable cut-off parameter value ~ 1.2 GeV. It turns out that the topological and 
collective quantum-wave fluctuations contribute to the QCD vacuum energy density with 
opposite signs. Therefore, a particular matching of numerical values of the non-perturbative 
QCD parameters (e.g. light hadron masses and scale parameters) via yet unknown dynamical 
mechanism could, in principle, provide zeroth net value of the non-perturbative QCD vacuum 
energy density without incorporating any extra physics at different space-time scales. 

Note, the hypothesis about the exact cancellation of quantum-topological and quantum- 
wave contributions to the vacuum energy as an internal feature of the theory of non- 
perturbative QCD vacuum does not mean that the sum of their quantum fluctuations is 



^ The numerical values for boundaries in Eq. (|2.3p can be somewhat mo del- dependent which, in practice, 
would not affect any of our conclusions here. 



24 




(2.4) 
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identically equal to zero. Indeed, let us consider the complete unordered two-point function 
(0|-i^".(x)Ff (x')|0) = (0| : (x') : |0) + (0|^i^",(x)Ff (x')|0)(.) . (2.5) 

TT Vr TT 



The first term in Eq. fl2.5p i s expressed via the experimentally measured value (introduced 
for the first time in Ref. |2^) and a correlation function D(x) which can be constrained e.g. 
in lattice QCD or effective field theory methods, i.e. j38| 

(0| : (a:') : |0) = (0| : ^i^",(0)Ff (0) : |0) Aop(x - x') , Aop(O) = 1 , 

vr vr 

(0| : — i^fc(0)F°(0) : |0) = (360 ± 20 MeV)^ . (2.6) 

TT 

For a detailed overview of the methods of calculation of higher power corrections to non- 



local condensates in QCD, see e.g. Ref. 39|. The second term in Eq. (12. 5p representing the 
quantum-wave component of the two-point function has an analogical representation. Now, 
the hypothesis about the exact cancellation of topological and quantum-wave contributions 
written in terms of normally ordered correlation functions 

(0| : ^FmK'W ■■ |0) = -(0|^i^»,(0)Ff (0)|0)(,) 

leads to the following expression for the sum of corresponding fluctuations: 

m-FUx)F:,\x'm = (0| : ^F^MK^O) ■. \0)D{x - x') , 

TT TT (2.7) 

D{x - x) = Dtop{x - x') - Dh{x - x') . D{0) = 0. 

As was mentioned earlier, the fluctuations of both types occur at the same space-time scale. 
However, identical vanishing of the function (12.71) can not be assumed here, i.e. a completely 
different space-time dynamics of the quantum-topological and quantum-wave fluctuations, 
of course, does not imply that D{x — x') = identically. 

Thus, we have illustrated another possibility to eliminate the non-perturbative QCD 
contribution to the ground state energy of the Universe based on the hypothesis about 
a specific fine-tuning of the QCD vacuum parameters. Both macroscopic previously dis- 
cussed in Ref. [29| and microscopic mechanisms considered above in this Section have a 
non-perturbative nature imprinted in essentially unknown quantum dynamics of the QCD 
vacuum, and the real situation can be, in principle, a superposition of both cancellation 
mechanisms (not excluding other possibilities, of course). 

Once the QCD vacuum energy contribution has been eliminated to the leading order 
in metric perturbations, the observable small A-term density can be further generated by 



quantum-gravity corrections according to the SZ scenario [19|, |20|], i.e. by quantum metric 
fluctuations dynamically induced by the non-zeroth non-perturbative gluon field fluctua- 
tions described by Eq. (12.71) . We will further demonstrate this fact by using the quasiclas- 
sical approximation methods in the conventional General Relativity and quantum gravity 
frameworks. 



III. QCD VACUUM ENERGY: FIRST ORDER IN METRIC PERTURBATIONS 
A. Equations of motion for metric fluctuations and macroscopic geometry 

In this Section, we briefly overview the quasiclassical (semiquantum) gravity framework 



in four dimensions, where the Zeldovich-Sakharov scenario [191 . |20[ | can be mathematically 
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realized in the simplest and well-grounded way. 

The semiciassicai approach to quantum gravity deals with quantum fields defined on a 
classical background |40|. Typically, one starts from the action of the gravitational and 
external physical fields written in terms of the corresponding quantum field operators as 
follows 

S = J Ld'x, L = -^^gg^^R,, + L{g'\xA) , (3.1) 

where g'^^ and Rik are the metric and curvature operators, respectively; L xa) is the 
Lagrangian density of external physical (e.g. gauge) fields xa with spins J < 2 interacting 
with each other and with gravity; index A numerates degrees of freedom of these external 
fields. 

In the framework of quasiclassical gravity theory it is conventionally assumed [40l] that 
the metric operator g^^ contains c-number part g'^^ - the macroscopic space-time metric, and 
operator part $f - the quantum graviton field. In what follows, we work in the Heisenberg 
representation. This means that one introduces an extra postulate about the existence of 
the Heisenberg state vector |0), which contains information about initial states of all incident 
quantum fields. Then, the graviton field operator by definition, satisfies the following 
condition on its average over the Heisenberg states: 

(0|$.» = 0. (3.2) 

As usual, derivation of the quasiclassical gravity theory equations of motion includes vari- 
ation and averaging operations. In order to match the classical evolution of the macroscopic 
metric with quantum dynamics of gravitons it is necessary to require that independent vari- 
ations of the action (13. ip over g^^ (taken at $f = const) and $f (taken at g"^^ = const) must 
lead to the same operator equations, namely, 

5 [ L(fx = -11 rf^x (^6g'^G,k) , =-l [d'^x (^5<^'^Gik) = , 

J ZJ\ /$*^=const Z J \ / c/»'==const 

(3.3) 

giving rise to 



El = ^ [g'm,^ - \5lg^'R,)j - 9'"%^ {g'\ xa) 



(3.4) 



where Tpm [g^^-, Xa) is the operator analog of the classical energy-momentum tensor corre- 
sponding to the Lagrangian L (^**'', xa)- By averaging the operator equations ( 13. 4p over the 
Heisenberg states, one obtains the equations of motion for the macroscopic (background) 
metric (7*^: 

(0|Gf|0)=0. (3.5) 

Subtracting c-number part (13.51) from operator (13. 4p . one obtains the equations of motion 
for the graviton fields : 

G,^-(0|G,^|0) = 0. (3.6) 

In addition to Eqs. (13. 5p and (13. 6p . one also gets the operator equations of motion for external 
fields with J < 2 by means of variations of the action (13. ip over xa- 
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An explicit form of the functional can be obtained by a variation of the action over 
the macroscopic metric g^^ without implying an explicit form for the quantum g^'^ {g''"^, 
operator. The condition (13. 3p providing consistency of equations (13. 5p and (13. 6 p fixes the 
exponential parameterisation for the metric operator as follows 41 



V^g'' = V^g'K^w = V^g'' (^f + + -^i^ri^L + ■■■), (3.7) 

where a shorthand notation for the graviton field has been introduced 



Applying the exponential parameterization (13. 7p . operator G\ can then be transformed to 
the following compact expression: 

+ ^ - l^'^R^ - , (3.8) 

where is the Ricci tensor of macroscopic space-time with the background metric gik in 
which all the covariant derivatives and lowering/ raising index operations are defined; and 

Tt = %) + \ ('^f '^r + g'"'gu) (^) ' g''f,m {g'\ xa) (3.9) 

is the total energy-momentum tensor operator incorporating the graviton field contribution 

1/1 



(3.10) 



Averaging the operator G\ given by Eq. (13. 8 p with an extra defining condition on the quan- 
tum graviton field $^ (13. 2p . we see that the equations of motion for the macroscopic (back- 
ground) metric (^jfc (13. 5p are transformed into usual Einstein equations as expected 

-(R\-\^\R\=Wm. (3.11) 



X V ' 2 



where a macroscopic average (0|Tj'^|0) in the right hand side contains a contribution from 
the energy-momentum tensor of the quantum graviton field T^^q-^ according to Eq. (13. 9p . 
Applying the same procedure to Eq. (13. 6p using expression (13. 8p . we finally obtain the 
equations of motion for graviton fields in explicit form: 

^ - V'!:? - ^ + + ^\R\ + ^\R\ - ^\RT^L = 2x (ff - (o|ff |o)) . (3.12) 

9 



A more rigorous derivation of the equations of motion (13.111) and (I3.12p is based on the 
canonical quantum gravity in the path integral (Faddeev- Popov) formulation j42|. In order to 
turn from the complete quantum gravity to its quasiclassical (semiquantum) limit, one starts 
from factorization of the path integral measure provided by the exponential parametrization 
of the metric operator (13.71) . A subsequent calculation of the factorized path integral, which 
is exact over the high-frequency (quantum) fields $^ and saddle-point approximated over 
the slow-changing (classical) fields Qik-, is equivalent to solving the equations of motion in 
the operator formulation (13.111) and (I3.12p . This demonstrates theoretical consistency of the 
operator approach described above. 



B. Operator gluodynamics with vacuum anomaly 

The non-perturbative fluctuations of gluon and quark fields naturally gravitate and should 
be included into equations of motion in the framework of operator field dynamics (13. lip and 
(I3.12p . The non-perturbative dynamics of gluon and quark fields is not sufficiently well 
developed in the literature, so we build up our analysis based upon the functional relations 
between metric fluctuations and vacuum fluctuations of quark and gluon fields only. These 
relations can be obtained from Eq. (13.121) and imply the existence of an adequate field- 
theoretical model for the QCD vacuum energy-momentum tensor incorporating conformal 
anomalies. The recipe for getting such a tensor is demonstrated e.g. in Ref. j35|: within 
the variational procedure for getting the energy-momentum tensor and equations of motion, 
the QCD coupling Qs can be viewed as an operator depending on operators of quantum 
fields according to the Renormalisation Group (RG) equations. In the framework of this 
formalism, one introduces the gluon field (vector-potential) operator A^, as a variational 
variable, related to the gluon field operator in the standard normalisation as follows = 
QsA-l. The stress tensor operator is then defined as J-]^ = diA^ — dkAi + f°'^'^A'lAl. In the 
case of pure gluodynamics, the QCD coupling operator g"^ = g1{J) depends upon the gauge 
operator of least dimension J = T'^^T^^ by means of the operator RG evolution equation 

^J^-^ = 9l{m9l{J)\. (3.13) 

where /9[5's( J)] is the QCD /3-function calculable in the standard quantum field theory frame- 
work. A solution of Eq. (I3.13P is the substituted into the effective operator Lagrangian 

whose variation w.r.t leads to the operator energy-momentum tensor of the gluon field 

-na -nkl \ /'Q 1 c;^ 

and operator equation of motion of pure gluodynamics 

^ah \ n A &\9l{J)\\ ^^k 
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After a straightforward covariant generalization, the operator gluodynamics defined by 
Eqs. fl3.15p and f l3.16p can be incorporated into the quasiclassical (semiquantum) gravity 
in the Heisenberg operator formulation given by Eqs. (13. lip and (I3.12p . In the framework of 
this formulation it is crucial that the energy-momentum tensor (I3.15P is conservative under 
the operator equations of motion (I3.16p . 
In the one-loop approximation we have 

where A is the QCD scale parameter discussed below; h = 6(0) = 11 is the one-loop (3- 
function coefficient of the pure gluodynamics (without quark fields). 

In order to construct a realistic operator energy-momentum tensor in QCD one has to 
incorporate quark fields. Formally, inclusion of the quark fields lead to extra terms in 
the energy-momentum tensor and to new operator equations of motion. However, this 
procedure can be simplified within the adopted phenomenological approach. It is clear from 
the beginning that an account for quark fields changes the numerical values of the /3-function 
coefficients. Also, non-perturbative quark-gluon fiuctuations happen at a characteristic scale 
of four-momentum transfers smaller than the double charm quark mass, so one has to fix 
b = 6(3) = 9 in Eq. fl3.17p and, correspondingly, in Eqs. f l3.15p and ( I3.16p . Further, an 
induced character of quark fiuctuations in the quark sea (effectively arising mostly due 
to gluon splittings in qq pairs at small momentum transfers) provides that QCD vacuum 
observables can be approximately expressed through the square of averaged gluon field 
fiuctuations. In particular, for quantum-topological quark-gluon fiuctuations it is well known 
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(0| : ss : |0) ~ (0| : uu : |0) = (0| : dd : |0) = -(0| : —F^^^F'^ : |0)Lg = -(225 ± 25 MeV)^ 

(3.18) 

where Lg ^ (1500 ± 300 MeV) ^ is the correlation length of fiuctuations which is normally 
calculated through the experimental data on quark and gluon condensates (and supported by 
the lattice QCD calculations). Its value is close to the minimal scale of fiuctuations given in 
Eq. (12. 3p at which their level is maximal. Within phenomenologically reasonable assumptions 
discussed above in Sect. [TTl the non-perturbative quantum- wave (hadron) fiuctuations occur 
at the same space-time scales, thus they should satisfy to a functional relation analogical 
to Eq. (I3.18p . Under these assumptions the operator relation between quark and gluon 
fiuctuations can be established by taking the trace of the averaged quark energy-momentum 
tensor and then applying its conservation condition and relations (12. 2 p and (I3.18P valid for 
both quantum-topological and quantum-wave contributions. This procedure provides us 
with the effective quark contribution to the QCD energy momentum tensor in the following 
form: 

T;(,),efr = ^(^« + ma + m,)fllg) . (3.19) 

Finally, adding Eqs. (I3.15P and (I3.19P written in the one-loop approximation gives phe- 
nomenologically motivated complete QCD energy-momentum tensor in operator form 

6eff = 6(3) + 8Lg{mu + + m^) ~ 9.6 , 
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and the operator equation of motion 



(^f In^) =0. (3.21) 

Of course, the resulting model expressions fl3.20p . f l3.2ip are approximate and restricted 
by quantitative phenomenological estimates (I2.2p and fl3.18p which characterize the non- 
perturbative quark-gluon fluctuations. 

Further simplifications are possible and make use of series expansion of the logarithmic 
operator function in small fluctuations as follows 

ej ^ e(o|j|o) J -mm 

(0|J|0) 

Note, the logarithmic function in the energy-momentum tensor fl3.20l) comes as a multiplier 
to an expression which has zeroth vacuum expectation value. Thereby, an account for the 
leading correction to the logarithm means that in the energy-momentum tensor, together 
with leading terms linear in 6J = J — (0| J|0), one also takes into account the higher-order 
terms in 6 J. Since the color factor suppresses the mean square fluctuation of the logarithm 
by a factor of 1/24, then replacement of the operator function under the logarithm by 
its averaged value can considered as a good approximation. Under such a replacement, 
g^ieJ) operator transforms into the usual QCD coupling at a characteristic scale of non- 
perturbative QCD fluctuations in the region of four-momentum transfers squared less than 
2 ^ (1500 MeVy, i.e. 



In — ~ 4 In 



Aqcd 

where Aqcd — 160 MeV is the QCD scale parameter. 

All subsequent calculations are controlled by the fact that the approximate expressions 
derived above must satisfy the energy-momentum tensor conservation. The replacement 
of the operator function under the logarithm by its averaged value transforms opera- 
tor equations fl3.2ip into the standard Yang-Mills equations. To this approximation, the 
energy-momentum fl3.20p should also transform (up to a constant multiplicative term and 
a constant additive term) into the standard Yang-Mills energy-momentum tensor. This 
result is achieved by a replacement of the conformal anomaly operator in the energy- 
momentum tensor fl3.20p by its averaged value. Then, turning back to original symbols 
A'} = Ai/gs, F°i^ = J^iklQs under the approximations adopted above, one writes finally 

^^(QCD) = ^ [-F^iF^' + In ^ - ^(Ol^F^FnO) . (3.22) 

^abpik ^ Q ^ j^ab^ ^abQ^ _ gjabc^c _ (3 23) 

As a matter of fact, we have not done anything radically new here - similar equations are 
often used in the Euclidean QCD framework where the vacuum energy-momentum tensor 
is estimated at instanton solutions of the classical Yang-Mills equations j32|. A proper 
covariant generalization of this framework will be applied for the A-term density calculation 
in the next Section. 
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IV. A-TERM CALCULATION 



Now, we have prepared everything what is needed for estimation of the dynamically 
induced A-term. In the right hand side of the macroscopic Einstein equations (13. lip the 
terms which correspond to the graviton- mediated interactions of the non-perturbative quark- 
gluon fluctuations are of the order of 0{G) and 0{asG) with G being the gravitational 
constant. Of course, in order to estimate the leading order effect, it makes sense to take 
into account only the first-order (linear) non-vanishing terms in gravitational constant G. 
Furthermore, due to an obvious smallness of the typical QCD space-time scales compared 
to the cosmological scales, the induced quantum fluctuations of metric should be considered 
at the Minkowski background. At last, it is sufficient to consider only the trace of the 
macroscopic Einstein equations R + 4x£:a = giving rise to the QCD-induced A-term 
density 

= -^(0|^ [f) fa'^^F^uFm + ^(0|T(a)|0) . (4.1) 

Let us start with expansion of the gluon stress tensor in series over the metric fluctuations 
(gravitons). From Eq. fl3.23p written in Riemann space 

(^"'^ - asr'-'Al^ v^g'^g'-'Fl = o , 

it follows immediately 

F-k = Ft, + \i^F^u - 4F^k - ^[F^i + 0{asG) , (4.2) 

where F^'^ is the usual stress tensor at the macroscopic background which does not account 
for interactions between gluon field and metric fluctuations (gravitons). The expansion 
(14. 2 p can be used only for extraction of the leading-order effect in the first term of Eq. (14. ip , 
which initially is of the order of 0{as). Thereby, higher terms 0{asG) in the expansion (14. 2 p 
generate corrections to the A-term density (14. ip of the order of 0(a^(j), which go beyond 
the one-loop approximation adopted here and therefore are omitted. 

The second term of Eq. (14. ip is formed by induced fluctuations of the metric (i.e. by a 
solution of equation (I3.12p ). which are unambiguously related with the gluon field fluctua- 
tions. To the leading order in gravitational constant G, in the left hand side of Eq. (I3.12p 
we keep only proper fluctuations of the quadratic form F^iF^^ without taking into account 
gravity. Then, using Eq. (I3.22p in the right hand side of Eq. (I3.12p . we obtain an important 
relation between the graviton field ipik and the gluon field strength F^^, to the respective 
order: 

- 4, f - <'/ + C;; = (-Ft^F^' + \sfF::^,Ff) in ^ . m 

TT V 4 J Aqcd 

At the next step one substitutes the expansion (14. 2 p into the first term of Eq. (14. ip . The 
second term in Eq. (14. ip can be calculated using the trace of the averaged energy-momentum 
tensor of gravitons (I3.10p (the averaging removes total derivatives due to symmetry of the 
Minkowski space-time background). After the averaging, the zeroth-order term in metric 
fluctuations (given by the unperturbed trace of the quark-gluon vacuum energy-momentum 
tensor) disappears due to the hypothesis about exact cancellation of quantum-topological 
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and quantum- wave contributions discussed in Sect. [TTl Implying Eq. (14. 3p . the resulting 
non-vanishing effect in the A-term density turns out to be quadratic in the graviton field 
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(0| ( v^^/^rV - 4^F„"^F^V^r) |0) 

V TT TT / 



^ :(0| Um, " - nt " - 2< " ) |0) • (4.4) 



16x^ ' 2 

The second term in Eq. (14. 4p can be identically transformed and simplified as follows: by 
transferring derivatives to one of the multipliers in each term one arrives at the differential 
form which can then expressed through the left hand side of Eq. (14.31) . Then applying the 
latter and combining all the terms together, the resulting A-term density takes a remarkably 
simple form convenient for phenomenological analysis: 

= -^In^^ (0|-^^^'' i^l - -M) |0) • (4-5) 

16 cAqcd TT \ A J 

The simplest way to obtain the physically meaningful result is to fix the Fock gauge 
ilj^.f^ = 0. In this case, according to Eq. (14.31) . ip = 0, so the final solution for the graviton 
field reads 

^P^ix) = x6eff In / d'x'Gix - x') ■ (^Ft,{x')F^\x') - 5',^F^,{x')FT\x')) . (4.6) 

Aqcd J ^ TT 47r / 

where Q{x — x') is the Green function satisfying the Green equation Q'\ = —6{x — x'). 
After substitution of Eq. (14. 6 p into Eq. (14. 5p . the respective averages should be calculated 
according to the following rules: 

{^\-n{^)F^\^) ■ (-FL{x')Fr{x') - {o\-FUx')Fr{x'm) m = 

TT \ TT TT / 

{0\^Ft,{x)FUxm{0\-F^\x)Frix'm + (4.7) 

TT TT 

{o\-FS{x)Frix'm{o\^F^\x)FUx'm, 

TT TT 

(0|-i^"(x)F,^„(x')|0) = '—{g.,g,m - mO(0| : -i^:,(0)F,"^(0) : \^)D{x - x') . (4.8) 

TT 96 TT 

The first rule, Eq. (14. 7p . demonstrates the exchange character of gravitational interactions 
of the gluon fluctuations, the first-order one in fluctuations (the zeroth-order single-point 
averages are explicitly subtracted). Also, Eq. (14. 8p . besides the gauge/Lorentz symmetry 
properties, incorporates the effect of compensation of the quantum-topological and quantum- 
wave contributions - the two-point function D{x — x') is defined as a difference between the 
corresponding correlation functions as reflected in Eq. (12. 7p . Substitution of these formulas 
into Eq. (14. 5 p leads to our final result for the QCD-induced dynamical A-term energy density: 



= --G(0| : ^i^^Ff : |0)^ x f ^) ' In ^ I d'yQ{y)D\y) 

TT V 8 y cAqcd Aqcd J 



(4.9) 



^1 ± 0.5) X lO^^^A MeV^ 
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Here, we introduced the dimensionless parameter 



A = 1 d'yg{y)D\y) , (4.10) 

which is defined in the Euchdean 4-space such that A > 0. Its numerical value has to be 
established by a dynamics in a complete non-perturbative QCD vacuum theory, or estimated 
in the static approximation e.g. in lattice QCD or in effective field theory approaches. 

In general, the expression (14. 9 p provides a good estimate for the dynamical contribution 
of the graviton-exchange interactions in the hadronic vacuum to the vacuum energy density 
of the Universe, it is unavoidable and should be taken into consideration in any theoretical 
or phenomenological analysis of the Dark Energy. 



V. DISCUSSION AND CONCLUSIONS 

The scenario for the A-term generation suggested above is realistic under two necessary 
and sufficient conditions: 

• the exact compensation of the non-perturbative quantum-topological and quantum- 
wave contributions to the QCD vacuum energy density in the zeroth order in quantum 
gravity effects; 

• the almost exact compensation of the graviton-exchange interactions of topological 
and wave fluctuations with a small residual effect summarized in Eq. (14. 9 p which can 
identifled with the observable A-term. 

The physical and mathematical premises for the flrst condition are contained in physically 
reasonable estimates for the energy density possibly coming from the flnite-time instanton 



solutions of the cosmo logical Yang-Mills flelds as discussed in Ref. [29[ and/or from vacuum 
hadronic wave modes given by Eq. (12. 4p . Both of these estimates show that the respective 
contributions compensating a large negative instanton energy (12. 2p . may have a similar 
nature being components of the same non-perturbative QCD vacuum. An interesting analogy 
of the co-existence of vacuum contributions of essentially different types, but of the same 
origin, takes place e.g. in solid state physics: in a crystallization process there are both 
negative (binding energy of atoms in the crystal lattice) and positive (zero-point collective 
fluctuations of the lattice itself) contributions to the "vacuum" energy of the medium. Of 
course, such an analogy is not exact and there is no exact compensation of these contributions 
in this case. 

The problem of exact compensation of vacuum contributions at the QCD energy scale 
seemingly has a dynamical nature. No doubts, the QCD vacuum state in the modern Uni- 
verse has been created during its cosmological evolution: the quantum-topological vacuum 
structures (e.g. instantons) have been created in real time around the quark-hadron phase 
transition in the Universe evolution. The quantum state of vacuum collective (wave) fluc- 
tuations of such structures should, therefore, be matched to the quantum state of these 
structures themselves. At the same time, the unprecedented accuracy in cancellation be- 
tween different vacuum contributions to many relevant digits (or their major flne tuning), 
if it takes place in Nature, forces us to consider such a matching of the vacuum topological 
structures and their collective fluctuations as a new physical phenomenon. Most probably, 
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an adequate description of this phenomenon is the subject of a new dynamical theory of the 
non-perturbative QCD vacuum which has not been yet created. 

The second condition above is tightly connected to the first one. Indeed, as is seen from 
Eqs. fl2.5p . ( 12. 6 p and f l2.7p . the representation of the two-point function D{x — x') as a differ- 
ence between topological and hadron correlation functions appears as a consequence of our 
assumption about the exact compensation of the topological and hadron wave contributions 
to the vacuum energy density calculated without the graviton-mediated corrections. As we 
will see below, the smallness of the QCD parameter A ~ 3 ■ 10~^ introduced in Eq. fl4.10p 
is required by phenomenological arguments - both by data on the A-term density itself and 
by QCD phenomenology. On the other hand, the parameter A is a functional of the differ- 
ence between the correlation functions defined at the same space-time scales (12. 3p . Thereby, 
there are physical and mathematical premises for the mutual cancellation of topological and 
wave contributions also in the effect graviton-mediated interactions, but such a cancellation 
is not exact. Thus, without a self-consistent theory of the non-perturbative QCD vacuum 
the actual value of the residual effect and, hence, the genuine value of A-term density can 
only be estimated by means of available experimental data and phenomenological arguments 
only. Based on purely phenomenological QCD arguments, let us try to move on towards a 
simple numerical estimation of the A-term density starting from Eq. (14. 9p . 

Indeed, the small parameter A ~ 3 ■ 10~^ can, in principle, be expressed through the well- 
known QCD parameters. Substituting only topological Dtop{x — x') or hadronic Dh{x — x') 
part of the complete correlation function D{x — x') into Eq. (I4.10p one immediately obtains 
A ~ 1. Meanwhile, it is known that naturally small QCD parameters are the light quark 
u,d,s masses, so let us assume that the value of (14.100 with the complete QCD correlation 
function D{x — x') is different from zero due to the chiral symmetry breaking effects. A small 
shift between the characteristic scales of topological and hadronic fluctuations induced by 
the chiral symmetry breaking can by given in terms of small current quark m^, m^, rUs 
masses 

|l/Ltop - 1/ Lh\ ruu + rud + rris , 

leading to a good estimate 

A = fc ■ ^ / T. > . (5-1) 

where ~ 1 is the dimensionless factor, and the factor l/(27r)'^ appears in the Fourier 
transform of the corresponding Green function. Then, the experimentally observable value 
of the A-term density (II. ip can be naturally obtained for k = 1 and the light quark masses 
satisfying the approximate sum rule + + ms — 100 MeV, which is in agreement with 
experimentally known values m„ = 1.5 — 5 MeV, rrid = 3 — 9 MeV, m<j = 60 — 170 MeV. 
In this case, the phenomenological formula for the A-term density, obtained based upon 
Eqs. (14. 9 p and (15. ip . contains all the basic QCD and chiral symmetry breaking parameters, 
i.e. 



a. 



TT 



-1 



.A = vrG(0| : : |0)^ X " , " M n In . (5.2) 



(27r)4 cAqcd a 



QCD 



In fact, Eq. (15. 2p provides a naive estimate for the dynamical QCD-induced positive con- 
tribution to the A-term energy density £a > which is close to the observed value within 
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a factor of few. For a better estimate, one would certainly like to evaluate the A-factor in 
Eq. (14. 9 P to a higher precision incorporating the chiral symmetry breaking effects e.g. in the 
lattice QCD framework or elsewhere. 

To conclude, we notice that under the above two conditions the A-term generation at the 
QCD energy scale appears to be a natural phenomenon for supersymmetric (or superstring) 
scenarios of high-energy physics (but would be confronted, perhaps, by Technicolor mod- 
els). Indeed, within these scenarios the vacuum subsystem of a quantum-topological nature 
exists only in the QCD sector. All other vacuum subsystems consist mainly of perturbative 
(weakly-coupled) vacuum fluctuations of fundamental fields and one may expect that e.g. 
a proper supersymmetry theory is capable of explaining the exact compensation of their 
contributions into the A-term density if J2hos ^hos = Sferm ^ferm satisfied. This can be 
seen e.g. from Eq. (II. ip where the respective energy scale (virtual particle mass of bosons 
Mbos and fermions Mferm) appears to the sixth power while bosons and fermions contribute 
to the vacuum energy with opposite signs. A detailed study of cancellations of this type is 
planned for a forthcoming study. 
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